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Objective: To evaluate the peripheral blood mononuclear cell (PBMC) expression levels 
of hemeoxygenase 1 (HMOX-V), superoxide dismutase 1 (SOD-1), vascular endothelial 
growth factor A (VEGF-A), transforming growth factor beta 1 (TGF-PV), interleukin 
(/L)-6, IL-15 and AdipoQ genes to study their association with preeclampsia (PE). 
Methods: A total of 177 pregnant women were recruited: 108 cases and 69 controls. 
Quantification of gene expression was measured by quantitative real-time polymerase 
chain reaction (PGR) using TaqMan probes. Results: Underexpression of VEGF-A and 
TGF-pi was a constant in most of the cases (80.91% and 76.36%, respectively) and their 
expression was associated with onset and/or severity of disease (p values < 0.05). IL-6, 
IL-15 and AdipoQ, showed low or no expression in PBMC samples evaluated. 
Conclusion: PBMC underexpression of VEGF-A and TGF-pi is a hallmark of PE in 
the study population. 

Keywords Preeclampsia, Gene expression. Mononuclear cells, Angiogenesis, Oxidative 
stress. 
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BACKGROUND 

Preeclampsia (PE) is included within a spectrum of conditions known as 
hypertensive disorders of pregnancy. PE is defined as hypertension and 
proteinuria and frequently occurs in the second trimester of pregnancy (1). 
Signs and symptoms of the disease include edema, headache, and in severe 
cases it is associated with damages to the liver and kidney, clotting 
abnormalities and intrauterine growth retardation (2). In Mexico, PE repre- 
sents the first cause of maternal death with a frequency of 25% among women 
in gestation stage and worldwide the frequency of the disease ranges between 
10% and 25% (3-5). 

The etiology of PE is still unknown, although genetic, immunologic and 
infectious risk factors have been described. The pathological pathways of 
disease involve an abnormal placental implantation process and consequently, 
oxidative stress and endothelial dysfunction (1,6—8). Due to the pathological 
processes involved, PE biomarker research has focused on the identification of 
the components (nucleic acids and proteins) implicated in these pathways. In 
this sense, placental and circulating mRNA/proteins of genes such as vascular 
endothelial growth factor A (VEGF-A), placental growth factor (PGF), 
transforming growth factor beta 1 (TGF-(31), hemeoxygenase 1 (HMOX-1), 
superoxide dismutase 1 (SOD-1), adiponectin (AdipoQ), interleukins 6 (IL-6) 
and 15 (IL-15), respectively, have been associated with PE, some with 
controversial results (9-16). Recent evidence shows that placental and 
trophoblast remains may be released into the peripheral blood of pregnant 
women and can be used for the identification of pathological pregnancies such 
as PE (17,18). Considering that these "peripheral artifacts" may reflect the 
pathological process of the preeclamptic placenta and consequently induce a 
change in gene expression in blood circulating cells, we evaluated the 
association between mRNA levels in human peripheral blood mononuclear 
cells (PBMCs) of seven genes (HMOX-1, SOD-1, VEGF-A, TGF-pi, IL-6, IL-15 
and AdipoQ) and PE in Mexican women. Additionally, their capacity to 
discriminate severity and onset PE was also tested. 

METHODS 

Patients and Samples Collection 

A total of 177 pregnant women with singleton pregnancy (108 PE cases and 69 
normotensive controls) were enrolled from the High Risk Pregnancy 
Outpatient Service of the Unidad Medica de Alta Especialidad (UMAE) # 23 
of the Instituto Mexicano del Seguro Social (IMSS) in Monterrey, Mexico, 
between September 2009 and April 2010. They provided written informed 
consent for their participation. Institutional Review Board approval was 
obtained (1905 Local Committee, ID number R-2010-1905-17) and the study 
was conducted in accordance with the Helsinki Declaration. 

The PE cases were diagnosed according to general guidelines of the 
International Society for the Study of Hypertension in Pregnancy (19). PE was 
considered as severe if the patient had a blood pressure > 160 mm Hg systolic 
or > 110 mm Hg diastolic on two occasions at least 6h apart while the patient 
was at bed rest, and a proteinuria of 5 g or more in a 24-h urine specimen, or 
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3+ or greater in two random urine samples collected at least 4h apart. The 
sub-classification of the cases group for early-late PE was made considering 
that additional maternal risk, such as maternal mortality is some 20-fold 
higher at less than 32 gestation weeks (GWs) than when PE occurs at term. 
Early PE was considered if present before 32 GW and late PE as that occurring 
with more than 32 GW (20,21). The control group consisted of normotensive 
women with uncomplicated pregnancies. Study exclusion criteria included 
chronic hypertension, chronic renal disease, gestational hypertension, gesta- 
tional diabetes mellitus, and underlying medical diseases such as rheumatoid 
arthritis and systemic lupus erythematous. All patients provided a peripheral 
blood sample at the time of diagnosis (PE cases) or when they entered the 
protocol (controls). The blood samples were collected using EDTA collection 
tubes and were processed quickly in a period not longer than 30 min. 

Gene Selection 

Selection of the study genes was made according to their association with PE 
and/or their central role on the pathologic processes of the disease cited in 
previous reports (7,9,10,12—18,22—25). Seven genes involved in oxidative stress 
(HMOX-1, SOD-1), angiogenesis-related events (VEGF-A, TGF-(31) and 
inflammation pathways (/L-6, 7L-15, AdipoQ), were included. 

RNA Isolation and cDNA Synthesis 

Mononuclear cells were purified by the standard Histopaque®-1077 (Sigma- 
Aldrich®, St. Louis, MO) protocol. Total RNA was isolated from cells according 
to the QiagenRNeasy® Mini Kit protocol (Qiagen, West Sussex, UK). RNA 
concentration and its quality were measured using a NanoDrop 
Spectrophotometer (NanoDrop Technologies, Wilmington, DE). cDNA was 
synthesized from 1 |.ig of total RNA using SuperScriptTM III First- Strand 
Synthesis Super Mix and random hexamers (Invitrogen, Carlsbad, CA) 
according to the manufacturer's indications. cDNA samples were stored at 
-20°C until use. 

Quantitative Real-Time PGR 

Validation of qPCR assays 

Efficiency, linear dynamic range of each quantitative real-time PGR (qPCR) 
assay, and the selection of the endogenous gene, were assessed by the 
construction of standard curves (26,27). In this stage, two pools of cDNAs were 
used: the first from normotensive women (controls); the second from women 
with a diagnosis of severe PE. cDNA from human microvascular endothelial 
cell line (HMEC-1) (28) was included as a positive control, in all the 
determinations (Supplementary file Si, shows the culture conditions for this 
cell line). Five serial dilutions of each cDNA were made in a concentration 
range from 345.4 ng to 2762.9 ng. All the qPCR reactions were performed in 
a final volume of 10 |.iL, IX of TaqMan Universal PGR Master mix and IX of 
the TaqMan Gene Expression Assay primer and probe mix. Each assay was 
made in triplicate for the target gene of interest {HMOX-1 — Hs01110250_ml, 
SOD-1 - Hs00916176_ml, VEGF-A - Hs99999070_ml, TGF-/31 - 
Hs99999918_ml, IL-6 - Hs00174131_ml, 7L-15 - Hs00174106_ml and 
AdipoQ — Hs006605917_ml) or for hypoxanthine phosphoribosyltransferase 
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(HPRTl), 18S ribosomal and glyceraldehyde-3-phosphate dehydrogenase 
(GAPDH) endogenous genes. qPCR reactions were performed on StepOne 
Plus Real-Time PGR System (Applied Biosystems, Foster City, CA) according 
to the manufacturer's instructions. 

To select the quantification method, the PGR efficiency of both, the gene of 
interest and each endogenous gene, was determined considering the variation 
of Gt values with template dilution (Gt gene — Gt endogenous) (AGt). In this 
step, Gt values were generated from an equivalent standard curve. 
Gonsequently, the AGt values were plotted against log input amount and the 
slope value was obtained from a semi-log regression line. The absolute value of 
the slope of AGT versus a log input <0.1 was considered as a criterion for 
passing the validation experiment. 

The choice of the endogenous gene was made considering three criteria: 
(1) The most consistent and uniform expression (Gt value) in the different 
cDNA pools used in the construction of standard curves; (2) amplification 
efficiency similar or equal to the interest gene according to the absolute value 
of its slope (slope of AGt versus log input <0.1); and relative abundance of its 
expression (more abundantly expressed than the target gene). 

Determination of gene expression level 

During the quantification of gene expression using the qPGR assay, 500 ng of 
cDNA from each patient sample were used in duplicate for the amplification 
round (according to the mixing conditions previously described). In each plate, 
a standard curve consisting of five points using cDNA from HMEC-1 cell line, 
in a concentration range of 32.32— 1034.2 ng was included in triplicate. Gt 
values were normalized to the HPRTl endogenous gene and gene expression 
was calculated as fold-change relative to the average of expression in the 
controls. 

Statistical Analysis 

Comparisons of risk factors and the clinical and personal characteristics 
among groups were performed using the Mann— Whitney (7-test for 
continuous variables and Chi-square or Fisher's exact test for categorical 
variables. Differences between numerical data of gene expression and PE 
severity or onset of PE were made using Student's t test or the Mann- 
Whitney ?7-test. The association between PE severity, PE onset and 
proportions of gene under/over expression was made using the Ghi-square 
test, considering a cutoff point of under- or overexpression of 20% in the 
expression rate of change (fold change >1.2) compared with controls. The 
odds ratio (OR) was calculated for significant comparisons. Statistical 
analysis was performed considering a significance level of 0.05 in Sigma 
Plot vll.O software. 

RESULTS 

A total of 177 pregnant women (108 PE cases: 60 mild and 48 severe PE, and 
69 controls) 15-40 years of age, with a singleton pregnancy between 21 and 
38 GW were included in the protocol. Thirty-seven percent of the cases were 
diagnosed with PE before 34 GW (early PE) and the remaining 63% had a late 
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Table 1. General characteristics and risks tactors for PE of the study population. 





Controis 


Miid PE 


Severe PE 




Characteristic 


(n = 69) 


(n = 60) 


(1 = 48) 


p Value 


Maternal age, median (range) 


26 (17-36) 


27 (15-40) 


27.5 (17-40) 


0.061 


Number of pregnancies, mean (range) 


2 (1-6) 


2 (1-9) 


2 (1-5) 


0.769 


Family history of PE, n {%] 


3 (4.35) 


3 (5.0) 


5 (10.42) 


0.364 


Family history of mellitus diabetes 


27 (39.13) 


27 (45.0) 


22 (45.83) 


0.713 


type 2, n (%) 










Primipaternity, n (%) 


21 (30.43) 


21 (35.0) 


12 (25.0) 


0.533 


Use of barrier contraceptive 


25 (36.23) 


19 (31.7) 


13 (27.1) 


0.578 


method, n (%) 










Nulliparous, n {%) 


25 (36.23) 


22 (36.7) 


19 (39.6) 


0.663 



Table 2. Endogenous gene selection. 



Endogenous gene 


Study group 


p Value 


Controis Ct mean (±SD) 


Severe PE Ct mean (±SD) 


18S Ribosomal RNA 


8.250 (±0.067) 


9.860 (±0.147) 


<0.001* 


GAPDH 


26.123 (±0.756) 


25.757 (±0.736) 


0.581 


HPRIl 


26.751 (±1.665) 


26.017 (±1.098) 


0.307 



*Stafistical difference. 



onset of disease (34-37 GW). Table 1 shows the general characteristics of the 
study groups. There were no differences in risk factors (such as use of barrier 
contraceptive method, family history of PE, maternal age or number of 
pregnancies) between cases and controls (p>0.05). 

The results of the qPCR assays validation, including the values related to 
endogenous gene selection, efficiencies of amplification and the choice of 
quantitative method, are summarized in Tables 2 and 3, respectively. GAPDH 
and HPRTl genes showed the most consistent and uniform expression (Ct 
value) in the samples used in the construction of standard curves (Table 2). 
GAPDH expression showed more stability (p = 0.581) than HPRTl (p = 0.307); 
both genes had Ct values closest to those obtained for the target genes, and 
were more abundantly expressed than them. Despite HPRTl slope values were 
lower than the GAPDH (Table 3), both genes are endogenous suitable for 
PBMC gene expression experiments. However, in this study, HPRTl gene was 
chosen as endogenous control because placental oxidative stress is a hallmark 
of PE (29,30) and GAPDH gene has a defined role in this process (31,32). 

Table 3 shows the absolute slope values obtained from semi-log regression of 
the ACT versus log input plot of VEGF-A, TGF-(31 and HMOX-1 genes, 
respectively (see Materials and methods section for details). Due to the fact 
that the slope absolute values were higher than 0.1 (criterion for passing the 
validation experiment), the quantitative method of the relative standard curve 
was selected over the 2"^^^* method for all genes evaluated in the study. 

According to the presence/absence of expression of the genes included in the 
study, the potential usefulness of HMEC-1 cell line for the construction of 
standard curves was evaluated. Our results showed that seven genes of 
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Table 3. Selection of the quantification metlnod. 



Gene Study group Endogenous gene ACt Slope absolute value 



VtGr-A Controls 


HPRIl 


2.69 


0.3067 




1 8S ribosomal 


Z 1 .40 


0. 1 466 




GAPDH 


4.89 


0.5065 


Severe PE 


HPRTl 


2.79 


0.1844 




1 8S ribosomal 


21.14 


0.3733 




GAPDH 


5.14 


0.4754 


TGr-pl Controls 


HPRJ] 


5.89 


0.2322 




1 8S ribosomal 


1 2.yo 


n 1 1 o > 
U. 1 1 o4 




GAPDH 


3.75 


0.3781 


Severe PE 


HPRT] 


6.02 


0.241 2 




IBS ribosomal 


12.38 


0.3170 




GAPDH 


3.79 


0.4091 


HMOX-1 Controls 


HPRT] 


1.30 


0.2151 




18S ribosomal 


17.52 


0.1013 




GAPDH 


0.84 


0.3610 


Severe PE 


HPRTl 


1.33 


0.2426 




18S ribosomal 


17.01 


0.3184 




GAPDH 


0.90 


0.4105 



interest were expressed in this cell line (Supplementary Table Si) and 
consequently its cDNA was included as a reference in the quantification stage. 

In our study, IL-6, IL-lb and AdipoQ genes, showed low (Ct values above 35) 
or no real time PGR amplification in the patients samples used in the 
validation assays. These results were consistent even in spite of having an 
increased quantity of cDNA (data not shown) and these genes were excluded 
from subsequent stages of the analysis. Figure 1 and Table 4 show the results 
obtained from evaluation of the expression levels of. HMOX-1, SOD-1, VEGF-A 
and TGF-pi genes. In general, all PE cases independently of severity and 
onset of the disease, showed uniform expression levels distribution of HMOX-1 
and SOD-1 genes, in a range of -4.1-2.6 for HMOX-1, and -3.6-5.5 for SOD-1, 
respectively; only an outlier expression value (—12.0) was observed for HMOX- 
1 gene (Figure 1). VEGF-A and TGF-pi underexpression was predominantly 
observed in an 80.91% and 76.36% of PE cases. Expression levels for these 
genes were ranged from —4.9 to 7.8 for VEGF-A gene, and from —37.6-35.5 for 
TGF-pi, respectively. The mean expression was 1.8 for VEGF-A and —10.2 for 
TGF-pi (Table 4). There were no differences in the expression levels of the 
study genes between PE subgroups (p>0.05). 

The association between PE severity and the expression of the interest 
genes was evaluated considering a rate of change of 20%. In this approxima- 
tion, the mild and severe PE cases were sub-classified according to their 
HMOX-1, SOD-1, VEGF-A and/or TGF-fil genes over/under expression (Table 
5). Patients with a rate of expression lower than 20% for these genes were 
excluded from the analysis. According to the data mentioned above 32.7% of 
the cases were excluded from HMOX-1 analysis, 31.8% were eliminated from 
SOD-1 evaluation and 10.3% of PE women were not considered for VEGF-A 
association tests. TGF-(31 expression rates were higher or equal to 20% in all 
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VEGF-A 



TGF-B1 



HMOX-1 



SOD-1 



Figure 1. Distribution of gene expression levels in the PE cases. Selected genes (VEGF-A, 
TGF-13], HMOX-1 and SOD-1, respectively) were analyzed by real-time PGR. Total PBMG 
RNA isolated from normal pregnant women (r) = 69) and vv'omen v^W'n PE (n= 108) were 
analyzed. Gene expression was calculated using thie relative standard curve mettiod 
and was expressed as told-chiange relative to ttie average of expression in ttie controls. 



the PE cases and therefore, all of them were included in the analysis. Table 5 
displays the results obtained from the association between PE severity and 
over/under expression level status oi HMOX-1, SOD-1, VEGF-A, and TGF-f31 
genes, respectively. Considering only cases with a minimal differential rate 
expression of change of 20% for the interest gene, the proportion of severe PE 
patients with overexpression of TGF-pi was higher (33.3%) than that observed 
for mild PE cases (15.3%) (OR=2.78; 95% CI: 1.1-7.0; p = 0.028). Expression 
levels oi HMOX-1, SOD-1 and VEGF-A genes were not associated to severity of 
the disease (Table 5). In a similar approximation, the association between PE 
onset and the gene expression level of the study genes was assessed (Table 5). 
Considering a rate of change of 20% only the overexpression of VEGF-A 
significantly increased the odds of early PE by 6.1 times (95% CI: 1.16-32.13) 
among the study population (p = 0.048). 



DISCUSSION 

PE is a pregnancy-specific hypertensive syndrome that causes substantial 
maternal and fetal morbidity and mortality (23). The origin of this pathology is 
unknown; however, both the theory on which clinical PE is a result of factors 
released by the placenta into the circulation, generating a widespread 
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Table 5. Analysis of PE cases with an expression fold cinange of 20% for the evaluated 
genes and their association with severity and/or onset of disease. 



Gene 


Expression status^ 




Severity (n = 


1 UO ) 




Onset (n = 


1 UO ) 


Mild PE Severe PE 


p Vaiuel 


Eariy PE Late PE 


p Value! 


HMOX-1 


Over 


12 


10 




6 


16 






Under 


26 


24 


0.842 


20 


30 


0.300 




Totalf 


38 


34 




26 


46 




SOD-] 


Over 


20 


22 




14 


28 






Under 


20 


1 1 


0.152 


15 


16 


0.194 




Total§ 


40 


33 




29 


44 




VEGF-A 


Over 


4 


5 




6 


2 






Under 


48 


39 


0.728 


29 


59 


0.048* 




Totali 


52 


44 




35 


61 




TGF-13] 


Over 


9 


16 




8 


17 






Under 


50 


32 


0.028* 


31 


51 


0.598 




Tofalft 


59 


48 




39 


68 





tFold change of 20%. 

fThirty-six, §thirty-five and ±twelve coses, respectively, did not show differences in their 
expression levels respect to controls and they were excluded from the analysis. 

tfOne patient was excluded from TGF-pi test because the cDNA sample was insufficient 
for qPCR amplification of this gene. 

Hp Value obtained by Chi square statistic. 

*Statistical difference. 

maternal vascular endothelial dysfunction, and the knowledge that some 
placental abnormalities and utero-placental ischemia may induce the detach- 
ment of placenta micro-particles into the mother's blood stream, converge in 
the common PE pathologic processes that include vascular injury, inflamma- 
tion and an increased oxidative stress. In this study, we test the possibility that 
placental factors or placental remains released into maternal circulation as a 
reflection of the pathological process of PE could induce a change in gene 
expression of seven genes in maternal PBMC from Mexican women. We found 
that at the mRNA level, maternal PBMC expression of 7L-6, 7L-15 and AdipoQ 
genes were not a useful tool to discriminate PE cases from normotensive 
controls because their quantification was close to the reliable limit for real time 
PCR detection. Compared with healthy pregnancies, IL-6, IL-15 and AdipoQ 
circulating proteins, have shown to be differentially expressed in complicated 
pregnancies including PE (14,33,34); however, their expression levels in PBMC 
from PE women have not been previously reported. These differences between 
studies are not unexpected and may be explained considering that protein 
quantification in previous reports has been determined in early stages of 
pregnancy when clinical PE was not present. The low PBMC expression of 
these genes in our study may be a consequence of pregnancy stages in which 
PE is already established; on the other hand, maternal circulating proteins 
may be a reflection of their release by other tissues including placenta and not 
a direct contribution of maternal PBMC. 

In this study, under- or overexpression of the oxidative stress-related genes 
HMOX-1 and SOD-1 in maternal PBMC were not a constant in PE cases, and 
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their usefulness to discriminate severity or PE onset was not observed. Some 
studies have reported an increased HMOX-1 and SOD-1 protein level in 
plasma of PE women related to controls. More recently, Nakamura et al., 
described lower mRNA concentrations of HMOX-1, SOD-1, HMOX-2 and 
catalase, in the cellular component of blood from Japanese PE patients and 
their mRNA values correlated with PE severity (13). Differences between this 
last report and our results may be related with differences between the studies 
designs and data analysis; for example, the blood cellular component used by 
Nakamura et al. in their protocol included all plasma-free blood components; 
it may include circulating trophoblast cells, or polymorphonuclear (such as 
neutrophils and eosinophils) cells. In our study, a PBMC separation using a 
density gradient was performed and therefore, our results represent only 
mRNA transcripts from blood cells within a specific density range. On the 
other hand, in this study, the expression level determination was made using 
the relative standard curve method whose data normalization involved an 
endogenous gene as a control of variations along sample manipulation. In the 
same sense, unlike plasmid constructs carrying the interest gene, we used 
HMEC-1 cDNA in the construction of standard curves and due to gene 
expression may be gestational-age dependent. A pool of RNA from normoten- 
sive controls (matched for GW) was used as a calibrator in data analysis. All 
the differences mentioned before could have an impact on the results obtained. 

In this study, the expression levels of both, VEGF-A and TGF-[31 genes were 
markedly diminished in the cases group; this finding was obtained without the 
sub-classification of cases by severity or onset of disease; these results agree 
with the central role of angiogenesis throughout pregnancy and its recognized 
imbalance in the development of PE. VEGF-A is the most ubiquitous member 
of the vascular endothelial growth factors family. It is a hypoxia-induced gene, 
and is positively regulated by TGF-pl and IL-6 in many cell types (10). 
Throughout normal pregnancy, vascular health is maintained through VEGF- 
A and TGF-pi binding to their respective membrane receptors (VEGF-Rl, 
VEGF-R2 and TGF-pRII, respectively) (11,35) on the endothelial cells. In PE, 
VEGF-A and TGF-pl pathways are antagonized by soluble proteins sFltl 
(soluble receptor of VEGF-A) and sEng (soluble co-receptor of TGF-pi). These 
anti-angiogenic soluble factors are a product of placental over secretion and 
cause local and systemic endothelial dysfunction (36). Many studies have 
reported increased plasmatic protein levels of sFltl, sEng, VEGF and TGF-pi, 
from pregnancies complicated by PE (10), but their expression levels in PBMC 
had not been described. Our results suggest that VEGF-A and TGF-pl 
production by PE PBMC may contribute in a small proportion in the increased 
protein level of these factors in maternal plasma, and it may be possible that 
elevated anti-angiogenic soluble proteins or other PE plasma components may 
have a negative effect on mRNA expression of VEGF-A and TGF-pl in PBMC 
from PE women. 

Unlike expression of HMOX-1 and SOD-1 genes, in this study, over- 
expression of VEGF-A was associated with an early onset of PE, and the 
overexpression of TGF-pi with severity of the disease. Even when we observed 
significant differences between the proportions of PE cases with overexpres- 
sion of VEGF-A or TGF-131 genes in the study subgroups, such proportions 
involved a small number of patients with this condition (8 PE cases with an 
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overexpression of VEGF-A and 25 with an overexpression of TGF-pi), 
therefore these association results should be considered cautiously. 

It is important to note that about 30% of PE cases showed a rate of change 
for HMOX-1 and SOD-1 very similar to that of the controls (lower than 20%), 
so it was not possible to include them in the association tests. Interestingly, 
only 10.28% of PE cases showed this feature for the VEGF-A, and it was not 
present in any patient for TGF-pi. We consider that this range of differential 
expression of angiogenic-related genes is an important finding that may 
ensure reproducibility of the results, and could be considered in gene 
expression approaches for PE biomarkers. 

In conclusion, our results show that the preeclamptic status is reflected on 
the expression level imbalance of VEGF-A and TGF-pi genes in maternal 
PBMC. The association between PBMC IL-6, IL-15 and AdipoQ expression and 
PE could not be assessed; therefore, their potential role in the pathology of 
disease must be evaluated in future approximations. 
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